Superpartner mass measurements with ID decomposed Mp 2 
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We propose a new model-independent technique for mass measurements in missing energy events 
at hadron colliders. We illustrate our method with the most challenging case of a short, single-step 
decay chain. We consider inclusive same-sign chargino pair production in supersymmetry, followed 
by leptonic decays to sneutrinos: x + X + ~ * t '■ ViV# (x~X~ ~* t~ t'~v^v\i). We introduce two 
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one-dimensional decompositions of the Cambridge Mt2 variable: Mt2^ and Mt2 ± , on the direction 

of the upstream transverse momentum Pt and the direction orthogonal to it, respectively. We show 
that the sneutrino mass M c can be measured directly by minimizing the number of events N(M C ) 
in which Mt2 exceeds a certain threshold, conveniently measured from the endpoint M^°^(M C ). 

PACS numbers: 14.80. Ly,12.60.Jv,11.80.Cr 



The Large Hadron Collider (LHC) at CERN has begun 
its long awaited exploration of the TeV scale, where new 
physics beyond the Standard Model (SM) may hold the 
key to our understanding of some very basic questions 
about our universe: What is the dark matter? What 
are the fundamental symmetries of Nature? Are there 
any hidden dimensions of space? A potential discovery 
of a missing energy signal at the LHC may relate to all 
three of these questions, if the missing energy is due to 
a stable, neutral, weakly interacting massive particle in 
a theory with space-time supersymmetry (SUSY) [l| or 
extra dimensions Q. 

The first order of business after the discovery of a miss- 
ing energy signal at the LHC will be to measure the mass 
of the missing particle and prove that it is not simply a 
SM neutrino [3( . This deceptively simple task turned out 
to be a notoriously difficult challenge. The generic topol- 
ogy of a prototypical "SUSY-like" missing energy event is 
schematically depicted in Fig. [1] Consider inclusive pro- 
duction of an identical pair of new particles P (from now 
on referred to as "parents"). Each parent decays semi- 
invisibly to a set of SM particles Vi, (i — 1,2), which 
are visible in the detector, and a dark matter particle C 
(from now on referred to as the "child" ) which escapes de- 
tection. In general, the parent pair may be accompanied 
by a number of additional "upstream" objects U (typ- 
ically jets) with total transverse momentum Pt. They 
may originate from various sources such as initial state 
radiation or decays of even heavier particles up the decay 
chain. We shall not be interested in the exact details of 
the physics responsible for U, adopting a fully inclusive 
approach to the production of the parents P. Given this 
general setup, the goal is to determine independently the 
mass M p of the parent and the mass M c of the child. 
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FIG. 1: The generic event topology under consideration. All 
particles visible in the detector are clustered into three groups: 
upstream objects U with total transverse momentum Pt, and 
two composite visible particles Vi, each with invariant mass 
mi and total transverse momentum pit- The transverse mo- 
menta of the two missing particles are labelled by p iT . 



p. In the simplest example of a short, single-step decay 
chain, each Vi consists of a single SM particle of fixed 
mass mi , and neither of these two approaches will work. 
One must then resort to methods based on the Cam- 
bridge Mt-i variable 0-i] or the related Sheffield M CT 
variable EH • Unfortunately, in order to apply those 
techniques, one must work with a subset of events within 
a relatively narrow fixed Pt range, incurring some loss 
in statistics. 

In this Letter we propose a new method which uses the 
full data set, with no such loss in statistics. Our method 
is based on the "subsystem" variant Q of the original 
Mt2 variable Q. For any given event, one can construct 
the transverse mass MiT of each parent P: 



Mi 
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where 



In the past, several approaches to this problem have 
been proposed, e.g. invariant mass endpoint measure- 
ments Q or exact reconstruction of the missing particle 
momenta p? T 5] . Unfortunately, they only apply to suf- 
ficiently long decay chains, where the visible particles in 
Vi arise from a sequence of at least three 2-body decays 
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is the transverse energy of the visible particle Vi and child 
particle C in each branch of Fig. [IJ correspondingly. The 
individual momenta pf T of the missing child particles C 
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are unknown, but they are constrained by the measured 
missing transverse momentum fx in the event: 



Pit + P2T 



-Pt 



Pit -P2T- 



(3) 



For the true values of the missing momenta p£ T , each 
transverse mass in ^ is bounded from above by the true 
parent mass Mp. This fact can be used in a rather in- 
genious way to define the Cambridge M?2 variable Q- 
One takes the larger of the two quantities in (fTJ) and 
minimizes it over all possible partitions of the unknown 
children momenta p£ T , subject to the constraint (j3|): 
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For a given Pt, the endpoint M^ aa: of this distribution 
gives the parent mass M p as a function of the input trial 
child mass M c : 



M p (M c ,P T ) 



M£T(M C ,P T ). 



(5) 



This property provides one relation among the two un- 
known masses M p and M c Q. 

Here we propose to obtain a second relation by using 
the property that the function M p (M c , Pt) is indepen- 
dent of Pt at the true child mass M c : 

M P (M C , P T + AP T ) - M P (M C , P T ) = 0, V AP T , (6) 

which we can rewrite more informatively as 

M P (M C ,P T )-M P (M c , 0) > , (7) 

with equality being achieved only for M c = M c . Eq. ([7]) 
implies that, for any given M c , there will always be a 
certain number of events whose Mt2 values will exceed 
the reference value M p (M c ,0), unless the trial mass M c 
happens to coincide with the true child mass M c . In 
order to quantify this effect, we define the function 



N(M C 



all events 
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where H(x) is the Heaviside step function. From the 
definition of N(M C ) it is clear that it is minimized at 
M c = M c , where in theory we would expect 



N„ 



i{N(M c )} = N(M c ) = Q. 



(9) 



In reality, the value of 2V m m will be lifted from 0, due 
to finite particle width effects, detector resolution, etc. 
Nevertheless we expect that the location of the N(M C ) 
minimum will still be at M c — M c , allowing a direct 
measurement of the child mass M c : 



{M c \N(M c )=N min }, (10) 



which is our first main result. Once the child mass M c is 
found from (|10[) . the true parent mass M p is obtained as 
usual from © as M p = M p (M c , P T ). 



At this point it is not clear whether we have gained 
anything statistics-wise, since the reference quantity 
M p (M c , 0) appearing in the definition ([5]) has to be mea- 
sured at a fixed Pt = anyway. Our second main result 
in this paper is that M p (M c , 0) can in fact be measured 
from the full data set with no loss in statistics as follows. 

Let us introduce one-dimensional (ID) decompositions 
of Mt2 onto the two special directions defined by the 
upstream momentum vector Pt- Following Ref. [ll|, first 
project the visible transverse momenta p^ of Fig. [T]onto 
the Pt direction (Tii) and its orthogonal direction (Tj_): 



PiT ± 



\s (St • Pt) Pt, (11) 
Pt x (p lT x P T ) , (12) 



PiT ~ PiT 



and similarly for the two transverse momenta p£ T of the 

children and for fr- Now consider the corresponding ID 
decompositions of the transverse parent masses |T]) 
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in terms of the ID projected analogues of ([2]) 
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Now we define ID Mt2 decompositions in complete anal- 
ogy with the standard Mt2 definition ((4]): 

M T2|| = min ^ {max {M iT|| , M 2 t ]] } } , (13) 
M T2± = min ^ {max{M 1Tx ,M 2Tx }}.(14) 

These decompositions are extremely useful. For once, 
the ID variables (|13I14[) can be calculated via simple an- 
alytic expressions as shown below. In contrast, a general 
formula for the original Mt 2 variable ([J) in the presence 
of arbitrary Pt is unknown and one still has to compute 
Mt2 numerically [HI]. More importantly, Mt 2± allows 
us to measure the reference quantity M p (M c ,0) in © 
from the full data set, using events with any value of Pr- 
To understand the basic idea, it is sufficient to consider 
the simplest, yet most challenging case of a single step 
decay chain. Let Vi be a single, (approximately) mass- 
less SM particle: mi = m 2 = 0. (The discussion for the 
massive case proceeds analogously.) In what follows, for 
illustration we shall use the same-sign dilepton channel in 
supersymmetry, where each Vi is a lepton resulting from 
a chargino decay to a sneutrino 0] . The charginos them- 
selves are produced indirectly in the decays of squarks 
and gluinos. For concreteness we shall use a SUSY spec- 
trum given by the LM6 CMS study point [13] . At point 
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FIG. 2: The unit-normalized Mt2 ± distribution (fl9j) for the 
same-sign dilepton channel in a SUSY model with LM6 CMS 
mass spectrum and a choice of test mass M c = 100 GeV. 
The yellow shaded distribution shows the theoretically pre- 
dicted shape (|19[1 . matching very well the parton level result 
from PYTHIA with no cuts (red histogram) . The green (blue) 
histogram is the corresponding result after PGS detector sim- 
ulation with mild (hard) cuts as explained in the text. The 
endpoint expected from eq. (|16|) is 132.1 GeV and is marked 
with the vertical arrow. 



LM6, the chargino (sneutrino) mass is M p — 305.3 GeV 
(M c = 275.7 GeV), and the rest of the SUSY mass spec- 
trum can be found in [l3| . In our simulations we use 
the PYTHIA event gen erator fu3 | and the PGS detector 
simulation program 

The variable Mt2± has several unique properties. 
Eventwise, it can be calculated analytically as 



T2 , 



An 



M} 



(15) 



A T ± = ^ ( \PlT ± I \P2T ± \+PlT ± ■ P2T ± ) • 

The endpoint of the Mt2±_ distribution is given by 



: (M C )= M +V M 2 + M 2 ) (16) 
in terms of the parameter \x introduced in Q 
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Eq. (|16[) reveals perhaps the most important feature of 
the Mt2 ± variable: its endpoint is independent of the 
upstream Pt and can thus be measured with the whole 
data sample. We can even predict analytically the shape 
of the (unit-normalized) differential Mt2± distribution 

^ r N 0± S(M T2± -M cW -N 0± )^-, (18) 

where Nq ± is the fraction of events in the lowest M c bin 



FIG. 3: Observable Mt2 ± distribution after hard cuts for 100 
fb" 1 of LHC data. The total stacked distribution consists of 
the SUSY signal (red) and the SM background (blue). The 
solid line is the result of a simple linear fit, revealing endpoints 
at 134.4 GeV and 172.4 GeV. 



Mt2± = M c , while the shape of the remaining (unit- 
normalized) Mt2 ± distribution is given by (see Fig. [5]) 
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Notice that this shape does not depend on any unknown 
kinematic parameters, such as the unknown center-of- 
mass energy or longitudinal momentum of the initial hard 
scattering. It is also insensitive to spin correlation ef- 
fects, whenever the upstream momentum results from 
production and/or decay processes involving scalar par- 
ticles (e.g. squarks) or vectorlike couplings (e.g. the QCD 
gauge coupling). It is even independent of the actual 
value of the upstream momentum Pt- Thus we are not 
restricted to a particular Pp range and can use the whole 
event sample in the Mt2 ± analysis. For any choice of 
M c (in Fig. H we used M c = 100 GeV), eq. dHJ) is a one- 
parameter curve which can be fitted to the data to obtain 
the parameter /i and from there the Mt2 ± endpoint (|16[) . 

As always, there are practical limitations to the use 
of such shape fitting. First, the shape (fl"§)) is modified 
in the presence of "mild" cuts, which are required for 
lcpton identification in PGS (green histogram in Fig. [2]), 
and more importantly, for the discovery of the same-sign 
dilepton SUSY signal over the SM backgrounds. To en- 
sure discovery, we use "hard" cuts as follows [HI, [l6| : ex- 
actly two isolated leptons with px > 10 GeV, at least 
three jets with p T > (175,130,55) GeV, P T > 200 
GeV and a veto on tau jets. With those cuts, in the 
dimuon channel alone, the remaining SM background 
cross-section is dominated by tt and is just 0.15 fb, while 
the SUSY signal is 14 fb, leading to a 22a discovery with 
just 10 ftT 1 of data [Hill. The distortion of the M T2 ± 
shape with these hard offline cuts is illustrated by the 
blue (rightmost) histogram in Fig. O The actual Mt2± 
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FIG. 4: The function N(M C ) denned in (22J). The blue (red) 
set of measurements are with (without) SUSY combinatorial 
background. The error bars shown are purely statistical. 

distribution which we expect to observe with 100 fb _1 of 
data, is shown in Fig. [3] and is comprised of a relatively 
small SM background component (blue) and a dominant 
SUSY signal component (red). In spite of the presence 
of a sizable SUSY combinatorial background, the Mt2 ± 
cndpoint expected from Fig. [3] is clearly visible and its 
location from a simple linear fit is obtained as 134.4 GeV, 
which is very close to the nominal value of 132.1 GeV. 
(Interestingly, the data reveals a second endpoint at 172.4 
GeV, which is due to events in which one chargino de- 
cays through a charged slepton: xf — >■ if, — > Xi \M- Its 
nominal value is 169.2 GeV.) 
Our final key observation is that 

M P (M C , 0) = M^ 2 ax (M c , 0) = M? 2 a *(M c ), (20) 
which allows to rewrite the function N(M C ) of eq. ([5]) as 
N(M C )= J2 h(m T2 -M^(Mc)). (21) 

all events 

The Mt2± analysis just described allows a very pre- 



cise measurement of the benchmark quantity M™ 2 ax (M c ) 
appearing in (|2"Tj) . so that the function N(M C ) itself can 
be reliably reconstructed, using the whole event sample 
all the way throughout the analysis, without any loss in 
statistics. We show our result in Fig. 21 where for conve- 
nience we unit-normalize the function N(M C ) as 

N(M C ) = N(M C )/(N(M C )), (22) 

where the averaging is performed over the plotted range 
of M c . As expected, the function N(M C ) exhibits a 
minimum in the vicinity of the true sneutrino mass 
M c = M c = 275.7 GeV. Ignoring the SUSY combina- 
torial background, this measurement (red data points) 
is quite precise, at the level of a few percent. In order 
to reduce the combinatorial background, we select events 
with M c < M T2± < M%?* and veto very hard[l7| lep- 
tons with pt > 60 GeV. The resulting M c measurement 
(blue data points) is at the level of 10%. This precision 
is clearly sufficient to exclude SM neutrinos as the source 
of the missing energy, hinting at a potential dark matter 
discovery at the LHC. 

In conclusion, we summarize the novel features and ad- 
vantages of our method in comparison to previous Mt2~ 
based proposals in the literature [1, 0. First, we make 
crucial use of property ([5]), which allowed us to measure 
directly the missing particle mass M c as in eq. (|10[) . Sec- 
ond, both the benchmark quantity M™ 2 ax (M c ) entering 

eq. (|2~Tj) as well as the the function N(M C ) itself can be 
measured using the whole available data sample at any 
Pt- To the extent that the definition of Mt2 ± relies only 
on the direction and not the magnitude of the upstream 
Pt, our method is insensitive to the jet energy scale error 
fill ] . We have also provided exact analytical formulas for 
the computation of the ID decomposed Mt2 variables [l8| 
and the shape (fT9"]) of the Mt2 ± distribution. 
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The measured value of M™^ in Fig. [3] already implies 
that the mass splitting M p — M c is on the order of 30 
GeV, resulting in a rather soft lepton pr spectrum. 
The corresponding analytical results for Mt2^ can be 
found in the first version of this paper, which is avail- 
able on the hep-ph archive. 



